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Charmonium Production in e^e Annihilation at = 10.6GeV 

Kuang-Ta Chao'' and Li-Kun Hao'^ 

'^Department of Physics, Peking University, Beijing 100871, People's Repubhc of China 

We review the charmonium production in e^e~ annihilation at ^/s — lO.GGeV. NRQCD predictions for the 
J/tp production are compared with the recent measurements by BaBar and Belle. A NRQCD calculation for the 
D-wave charmonium production is reported to further test the color-octet mechanism and detect the J^'^' = 2 
state. The issue of double cc production in e~^e~ annihilation observed by Belle is also discussed. 



1. INTRODUCTION 

The newly developed nonrelativistic QCD 
(NRQCD) factorization formalism ^ allows 
the infrared safe calculation of inclusive heavy 
quarkonium production and decay rates. In the 
NRQCD production mechanism, a heavy quark- 
antiquark pair can be produced at short distances 
in a conventional color-singlet or a color-octet 
state, and then evolves into an observed quarko- 
nium nonperturbativcly. With this color-octet 
mechanism, one might explain the Tevatron data 
on the surplus production of J/ ip and ip' at large 
Pt, but puzzles about their polarizatons still re- 
main (for a review see Q and references therein) . 

To further test the color octet mechanism, it 
may be interesting to study the charmonium pro- 
duction in e~^e~ annihilation. The J/ip produc- 
tion in e'^e~ annihilation has been investigated 
within the color-sing let model |,|,| and the 
color-octet model (6|j^,^. The angular distribu- 
tion and energy distribution of color-singlet J/tp 
production at ^/s = lO.QGeV have been discussed 
in ||]. In 1^ it is noted that a clean signature of 
the color-octet mechanism may be observed in the 
angular distribution of J/iJj production near the 
end point region. In |Q contributions of the color- 
octet as well as color-singlet to the J/ip produc- 
tion cross sections are calculated in a wide range 
of e~''e~ collider energies and it is found that with 
reasonable choices of both color-singlet and octet 
matrix elements the color- octet contribution will 
dominate. Moreover, the J/ip polarizations are 
predicted in Recently, BaBar ^ and Belle 
|lO[ have measured the direct J/ip production in 



contimmm e^e^ annihilations at y/s ~ lO.BGel/. 
The total cross section and the angular distribu- 
tion seem to favor the NRQCD calculation over 
the color-singlet model but some issues still 
remain. The P-wave charmonium XcJ production 
in e+e~ annihilation has been discussed in [ |l3| . 
The total Xc\,i cross sections are predicted to be 
dominated by the color-octet process. In addi- 
tion, as a further test of the production mecha- 
nism, we will report a calculation for the -D-wave 
charmonium production in ef^ e~ annihilation, to 
which the color-octet process will make substan- 
tial contributions Q, and these states may be 
detected in the future by BaBar and Belle. 
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The leading order color-singlet contributions to 
direct J/ip production include the following pro- 
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Eq.(|l|) is the gluon process, and Eq.(g) is the 
quark process (the charm quark fragmentation), 
and they are shown respectively in Fig. 1(a) and 
Fig. 1(b). Eq.(||) is the gluon jet process with 
q ^ u, d, c quarks. 

At y/s ~ 10.6GeV the gluon process (1) dom- 
inates over the other two color singlet processes. 
The process (3) is negligible at low energies, but 
grows as the energy increases, and at high enough 
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Figure 1. The main Feynman diagrams for the 
charmonium production in e^e^ annihilation, 
(a) and (b) are color-singlet processes; while (c) 
and (d) are color-octet processes. 



Figure 2. Color-singlet energy distribution da / dz 
as a function of z = —fi- at Jl — \0.&GeV. 
The distributions are from the giuon process 
(1) (dashed line), the quark process (2) (dotted 
line) , along with the sum of them (solid line) . See 
Ref.0. 

energies {^/s > 250GeV) its contribution will 
dominate over processes (1) and (2) [Q. 

In Fig. 2, We display the energy distributions of 
the gluon process (1), the quark process (2), and 
their sum at = 10.6 GeV. The energy spec- 
trum of these two processes are both flat and the 
gluon process contribution is much larger than 
the quark process (see 0). 

The leading order color-octet contributions to 
direct J/ip production in e~^e~ annihilation in- 
clude the following two processes: 

e+e- ^ r ^9 + cm^'+'Lj], (4) 
e+e- ^ 7* ^gg + cc[8,3 5i], (5) 

as shown respectively in Fig. 1(c) and Fig. 1(d). 
Here '^^'^^Lj denotes the states ^Sq and ^Pj, q 
represents the u, d, s, c quarks. At low energies 
the dominant process is the process (4), whereas 
the process (5) will become dominant at y^s > 
25GeV. 

As shown above, with a reasonable range of 
the color-singlet and octet matrix elements, for 
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the color-singlet cross section the quark process 
(2) contributes 0.07-0.1 pb §, which is consistent 
with the gluon process (1) contributes 0.2-0.3 
pb, and the total color-singlet cross section is 0.3- 
0.4 pb. The color-octet cross section contributes 
0.6-0.7 pb. Including both the color-singlet and 
octet contributions, we get a total cross section 
of about 0.9-1.1 pb. BaBar gives the total cross 
section cr(e+e" J/^jX) = 2.52±0.21±0.21 p6, 
and Belle gives cr(e+e" J/^X) = 1.47±0.10± 
0.13 pb. Their values are much larger than the 
predicted color-singlet value of 0.3-0.4 pb. More- 
over, for the J/V' angular distribution parame- 
ter A, at high momentum p* , NRQCD predicts 
0.6 < A < 1.0 1^ while the color-singlet model 
predicts A k, —0.8, BaBar gives A = 1.5 ± 0.6 for 
p* > 3.5GeV, clearly favoring NRQCD. 

However, there are some problems for the color- 
octet model. First, the expected peak at the up- 
per end point in the J/^ momentum p* distribu- 
tion in process (4) has not been observed. It is not 
clear whether this could be due to the smearing 
effect of soft gluon emission at the end point [jll] . 
Second, the J/ijj polarization observed by BaBar 
and Belle is mainly longitudinal, but the process 
(4) predicts the unpolarized J/t/i [1). The color- 
singlet process (1) predicts longitudinally polar- 
ized J/ip but it only gives a small portion of the 
cross section (~ 0.3 pb). Again, the polariza- 
tion problems remain as at the Tevatron for the 
NRQCD color-octet mechanism One of the 
possible solutions for the J/ijj polarization puz- 
zle is that the color Ml transitions, which will 
flip the spin of the charm quark and give rise to 
longitudinal polarization, may not be unimpor- 
tant as compared with the El transitions in the 
evolution of the charm quark pair. 

3. Xcj PRODUCTION IN e+e ANNIHI- 
LATION 

The P-wave charmonium XcJ production in 
e+e" annihilation has been discussed in [^. Be- 
cause the C-parity forbids the process e+e~ 
cc['^Pj]gg, the color-singlet processes are e+e~ 
XcJ + cc and e"''e~ — > XcJ + 7; a-nd the color- 
octet processes are e+e^ {cc)s[^ Sq,^ Pj] + 
g. We get these cross sections as a{e'^e^ 



Xc.n) = 0.001, 0.012, 0.005 p6, a{e+e- ^ Xcj + 
cc) = 0.005,0.018,0.008 and aoctetiXcj) = 
0.022,0.067,0.112 pb for J = 0,1,2 respectively. 
Then we get their total cross sections 

a{e+e- XcjX) = 0.028, 0.097, 0.125 ph, (6) 

which are consistent with |l^. The color-octet 
contributions are larger than the color-singlet by 
a factor of 2 ~ 10. These total cross sections are 
smaller than the upper limit given by Belle. 

4. D-WAVE CHARMONIA 5,, PRODUC- 
TION IN e+e ANNIHILATION 

In a recent study we find that the spin 
triplet D-wave charmonia 5j{J ~ 1,2,3) pro- 
duction in e+e~ annihilation at ^/s — 10.6 GeV 
mainly comes from Fig. 1(a) for the color-singlet 
process, and Fig. 1(c) for the color-octet pro- 
cess. The contribution of color singlet process 
in Fig. 1(b) to the 5i cross section is found to be 
smaller by a factor of about 60 than that for the 
J /ip and therefore is negligible, which is consis- 
tent with . Using the velocity scaling rule and 
relating the 5j color-octet matrix elements to that 
of t/"', we get the total cross sections 

a{5i) = CTi erg - 0.027 -I- 0.016 = 0.043p6 
a{52) = cTi + fjg - 0.067 -f 0.027 = 0.094p6 (7) 

However, if we use a more radical choice of the 
5 J color octet matrix elements by relating them to 
that of J/'f/', we would get much larger values for 
the 5 J production cross sections, cr((5i) — 0.16p6 
and cr((52) — 0.29p6. Therefore the measurement 
of 5j production in the future will be helpful 
to clarify the color-octet mechanism and the 5j 
color-octet matrix elements. 

In particular, for the 2 D-wave charmonium 
82 , which is expected to be below the open charm 
DD* threshold and forbidden to decay to DD, 
its width is estimated to be 300 — AOQKeV and 
its branching fraction of decay mode J/i/jtt+tt" 
is B{52 J/ijTT+TT-) ~ 0.12j|6), which is only 
smaller than that of B{ip' — > J/-07r+7r~) = 
0.324 ± 0.026 by a factor of 3. With more data 
available at the B Factories in the future, it will 
be possible to detect this 2^^ D-wave charmo- 
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nium state. The 1 D-wave state ip" (3770) 
could also be detected via ip" DD decay. 



We also thank K.Y. Liu, C.F. Qiao, and F. Yuan 
for their collaboration. 



5. PUZZLE OF DOUBLE cc PRODUC- 
TION IN e+e ANNIHILATION 

Recently the Belle Collaboration has reported 
a measurement on the double cc production in 
e+e" annihilation at ^/s = 10.6GeV, and found 
that a very large fraction of the produced J/ ip is 
due to the double cc production in e+e^ annihi- 
lation 



a{e 



J/ipcc)/a{e+e- J/ipX) 
= 0.59t°:i5±0.12. 



(8) 



which corresponds to cr(e+e~ —^ J/ipcc) ~ 0.9pb. 

This result is puzzling in terms of perturba- 
tive QCD calculations. The color-singlet process 
(2) (the charm quark fragmentation, see Fig. 1(b)) 
will contribute only ~ 0.1 pb to the cross section. 
Moreover, the produced J/ip should be trans- 
versely polarized in process (2)§, in contrast to 
the observed J/ijj being mainly longitudinally po- 
larized [p|JlOt|. The color-octet process (5) may 
also contribute to the double cc production with 
q = c, but again its cross section is too small. 
Whether the QCD radiative correction to process 
(2) can enhance the double cc production by a fac- 
tor of ~ 9 remains interesting, but it may still be 
difficult to explain the observed J/ ip polarization. 
So, we intend to conclude that it is very hard to 
explain the double cc production data observed 
by Belle based on NRQCD or, more generally, on 
perturbative QCD, and possible nonperturbative 
QCD effects have to be considered. 

In summary, we find that charmonium produc- 
tion in e"'"e~ annihilation at ^/s = 10.6GeV is 
interesting in testing the heavy quarkonium pro- 
duction mechanism in NRQCD. We expect that 
with more data collected at BaBar and Belle in 
the future, many theoretical issues will be further 
clarified, and we will reach a new step towards 
a better understanding for heavy quarkonium 
physics in QCD. 
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